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Matrix element UtN of sterile neutrino N mixing with Vt is the least constrained in the lit- 
erature among the three UaN [a = e, r) mixing parameters characterizing the sterile neutrino 
phenomenology. We study the contribution of massive dominantly sterile neutrinos to purely lep- 
tonic T-decays and semileptonic decays of t and K, D mesons. We consider some decays allowed 
in the Standard Model (SM) as well as Lepton Flavor and Lopton Number Violating (LFV, LNV) 
decays forbidden in the SM. From the existing experimental data on the branching ratios of these 
processes we derived new limits on Utn more stringent than the ones existing in the literature. 
These limits are extracted in a model independent way without any ad hoc assumptions on the 
relative size of the three different sterile neutrino mixing parameters. 
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I. INTRODUCTION 

Lepton flavors are conserved in the Standard Model (SM) due to the presence of an accidental lepton flavor 
symmetry, which, however, is broken by non-zero neutrino masses. Neutrino oscillation experiments have proven 
that neutrinos are massive, although very light, particles mixing with each other. Moreover, neutrino oscillations is 
the first and so far the only observed phenomenon of lepton flavor violation (LFV). In the sector of charged leptons 
LVF is strongly suppressed by the smallness of neutrino square mass differences {rnl. — mf,.)/qQ compared to the 
characteristic momentum scale, Qq, of an LFV process which is typically of the order of the charged lepton mass 
qq ^ mi. If neutrinos are Majorana particles there can also occur lepton number violating (LNV) processes. They 
are also suppressed by the smallness of the absolute value of m^. However, the situation may dramatically change if 
there exist either heavy neutrinos Ni, known as sterile, mixed with the active flavors Ve,^i,r or if there are some new 
LFV and LNV interactions beyond the SM. 

Here we study the former possibility and consider an extension of the SM with right-handed neutrinos. In the case 
of n species of the SM singlet right-handed neutrinos v'pjj = [v'fn, besides the three left-handed weak doublet 

neutrinos u'j^^ = {v'j^^, '^l^^ ^Lt) neutrino mass term can be written as 
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= -2C^'m^yiVi + ^m^M''^Vj) + h.c. (2) 
i=i i=i 

Here M-Lt-M-r are 3x3 and nxn symmetric Majorana mass matrices, and Md is a 3 x ti Dirac type matrix. Rotating 
the neutrino mass matrix to the diagonal form by a unitary transformation 

C/^X('')C/ = Diag{m,i,---,m,3^„} (3) 

one ends up with 3 -I- n Majorana neutrinos with masses m^^^ , ■ • • , ^1,3^^ . The matrix Uak is a neutrino mixing matrix. 
In special cases among neutrino mass eigenstates there may appear pairs with masses degenerate in absolute values. 
Each of these pairs can be collected into a Dirac neutrino field. This situation corresponds to conservation of certain 
lepton numbers assigned to these Dirac fields. Generically in this setup neutrino mass eigenstates can be of any mass. 
For consistency with neutrino phenomenology (for recent review, cf. [T]) among them there must be the three very 
light neutrinos with different masses and dominated by the active flavors Va (ck — e, /i, t). The remaining states may 
also have certain admixture of the active flavors and, therefore, participate in charged and neutral current interactions 
of the SM contributing to LNV and LFV processes. Explanation of the presence in the neutrino spectrum of the three 
very light neutrinos requires additional physically motivated assumptions on the structure of the mass matrix in ([I]) . 
The celebrated "see-saw" mechanism [2j , presently called type-I see-saw, is implemented in this framework assuming 
that ^ A^£)- Then, there naturally appear light neutrinos with masses of the order of ^ dominated 
by Va- Also, there must be present heavy Majorana neutrinos with masses at the scale of ^ Mr- Their mixing with 
active neutrino flavors is suppressed by a factor ~ AAd/M-r which should be very small. In particular scenarios 
this generic limitation of the see-saw mechanism can be relaxed [3]. Then the heavy neutrinos could be, in principle, 
observable at LHC, if their masses are within the kinematical reach the corresponding experiments. Very heavy or 
moderately heavy Majorana entry Mr of the the neutrino mass matrix naturally appears in various extensions of the 
SM. The well known examples are given by the S'O(10)-based supersymmetric jl] and ordinary [S] grand unification 
models. The supersymmetric versions of see-saw are also widely discussed in the literature (see, for instance, [6] and 
references therein). 

In the present paper we study the above mentioned generic case of the neutrino mass matrix in ([T]) without implying 
a specific scenario of neutrino mass generation. We assume there is at least one moderately heavy neutrino N in 
the MeV-GeV domain or even lighter. The presence or absence of these neutrino states, conventionally called sterile 
neutrinos, is a question for experimental searches. If exist, they may contribute to some LNV and LFV processes as 
intermediate nearly on-mass-shell states. This would lead to resonant enhancement of their contributions to these 
processes. As a result, it may become possible to either observe the LNV, LFV processes or set stringent limits on 
sterile neutrino mass m^r and mixing UaN with active neutrino flavors Va (ct — e,/i, t) from non-observation of the 
corresponding processes. 

On the other hand the sterile neutrinos in this mass range are motivated by various phenomenological models [7], 
in particular, by the recently proposed electroweak scale see-saw models [8], [9]. They may also play an important 
astrophysical and cosmological role. The sterile neutrinos in this mass range may have an impact on Big Bang 
nucleosynthesis, large scale structure formation (TU], supernovae explosions [TT]. Moreover, the keV-GeV sterile 
neutrinos are good dark matter candidates p!^HH| and offer a plausible explanation of baryogenesis [TS] . Dark Matter 
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sterile neutrinos, having small admixture of active flavors, may suffer radiative decays and contribute to the diffuse 
extragalactic radiation and x-rays from galactic clusters [16]. This is, of course, an incomplete list of cosmological 
and astrophysical implications of sterile neutrinos. More details on this subject can be found in Refs. [T7], [15] . 

The phenomenology of sterile neutrinos in the processes, which can be searched for in laboratory experiments have 
been studied in the literature in different contexts and from complementary points of view (for earlier studies see [T^ ) . 
Their resonant contributions to r and meson decays have been studied in Refs. [20l - f25] . Another potential process to 
look for sterile Majorana neutrinos is like-sign dilepton production in hadron collisions |26H29) . Possible implications 
of sterile neutrinos have been also studied in LFV niuonium decay and high-energy muon-electron scattering [30j . 
An interesting explanation of anomalous excess of events observed in the LSND [3T] and MiniBooNE [32] neutrino 
experiments has been recently proposed [33] in terms of sterile neutrinos with masses from 40 MeV to 80 MeV. An 
explanation comes out of their possible production in neutral current interactions of i/^ and subsequent radiative 
decay to light neutrinos. 

Here we study a scenario with only one sterile neutrino state N. Phenomenology of a single sterile neutrino N is 
specified by its mass niN and three mixing matrix elements UeN, U^^n, UtN- In the present paper we focus on the 
derivation of limits on the matrix element UTN^ which is currently least constrained in the literature. Towards this 
end we use the results of experimental measurements of branching ratios of purely leptonic r decays and semileptonic 
decays of t and K, D mesons |34j . One of the key points of our derivation is its model independent character, in the 
sense that we do not apply any additional assumptions on the relative size of the three mixing parameters UaN- Such 
ad hoc assumptions are typical in the literature and stem from the fact that all these three parameters enter in the 
decay rate formulas of any decay, potentially receiving contribution from N as an intermediate state. Therefore, in 
order to extract individual limits on each mixing parameter one may need additional information on them. We will 
show that in purely leptonic t decays it is unnecessary and in the other cases this sort of information can be procured 
by a joint analysis of certain sets of leptonic and semileptonic decays of r and K, D. 

The paper is organized as follows. In the next Section [ll] we present decay rate formulas for t and pseudoscalar 
meson LFV and LNV decays in the resonant domains of sterile neutrino mass rnjv- In Section III we derive upper 
limits on IC/twI from the existing experimental data on purely leptonic 5-body r decays, semileptonic r and K,D 
decays, considering sterile neutrino contribution as an intermediate state and in some cases as one of the final state 
particles. Section IV contains summary and discussion of our main results. 



II. DECAY RATES 



Neutrino interactions are represented by the SM Charged (CC) and Neutral Current (NC) Lagrangian terms. In 
the mass eigenstate basis they read 

= 4 E f^'^ + E ^"^f^- ^^^rPLl^3 Z^, (4) 

where I = e, and i = 1, ...,n + 3. We consider the case with a single sterile neutrino N and, therefore, we choose 
n = 1 and identify N — v^. 

In what follows we study sterile neutrino contribution to the following decays 

T- rere+vv, T- ^r^-K^TT- , M+^l+l^TT^ (5) 
n^N, -> rPiN, (6) 

where M — K, D, B and I, li — e, /i. In the first decay of Eq. ([S]) both v denote the standard neutrino or antineutrino 
dominated by any of the neutrino flavors Ve,Vp,,i'T- These reactions include lepton number and flavor conserving as 
well as LFV and LNV decays. In the first case they receive the SM contributions, which alone give good agreement 
with the experimental data. 

The LFV and LNV decays (l5| are only possible beyond the SM. In the present framework they proceed according 
to the diagrams shown in FiglT] with sterile neutrino as a virtual particle. Considering LNV decays we assume 
that sterile neutrino is a Majorana particle N — N'^. When the intermediate sterile neutrino A^ in these diagrams is 
off-shell their contribution to the processes ^ is negligibly small [2T], being far away from experimental reach. On 
the other hand there exist specific domains of sterile neutrino mass mpf where N comes, for kinematical reasons, close 
to its mass-shell leading to resonant enhancement [2DH^ of the diagrams in FigjI] These domains of will be 
specified below. 

The decay rate formulas for the reactions in Eq. ([s]) can be directly calculated from the diagrams in Fig [l]and 
Lagrangian Q for arbitrary mass mjv of sterile neutrino. We focus on the regions of m^v where the sterile neutrino 
contribution is resonantly enhanced |20H22j . In these mass domains the intermediate sterile neutrino in Fig. fl] can 
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FIG. 1; Structure of the lowest order contribution of sterile neutrino N to the semileptonic and leptonic meson and r decays. 
Here U, Ij — e, fi. 



be treated as nearly on-mass-shell state. This is to say, the sterile neutrino N is produced in the left vertices of 
the diagrams in Figjl] propagate as a free unstable particle and then finally decays in the right vertices. Thus the 
decay rate formulas for the reactions r, M — >■ X1X2I can be represented in the form of products of the two factors: 
T or meson decay rate to the sterile neutrino r(T, M — > NXi) and a branching ratio of the sterile neutrino decay 
Br{N — ^ 1X2) ^ where Xi,l represent final state particles of ([5|. This representation is approximate and valid in the 
"narrow width approximation" Fjv <SC where Fjv is the total decay width of sterile neutrino. As seen from Fig|6j 
this condition is satisfied in the region oi tum studied in our analysis where F^f < 10"^" MeV . Below we list the 
decay rate formulas in this approximation for the reactions in Eq. ([s]) specifying the corresponding resonant regions 
oi rriN where they are applicable. These formulas are readily derived from the diagrams in Figjl] considering the two 
vertices as the two independent processes of sterile neutrino production and its subsequent decay. 
For semileptonic decays of mesons M and r-lepton the decay rate formulas are 



F(M+ ^ 7r"e+e+) 
F(M+ ^ 7r-^+e+) 
F(M+ ^ ^l+^^e+) 
valid in rUe + < < itim — rrie, 

V{t- TT-Tl^V^) W F(t 



F(M+ ^ l+N) 
F(M+ ^ e+N) 
T{M+ e+N) 



F(iV'= ^ e+TT-) 

TV ' 
F(iV'= ^ M+TT-) 

TV 

F(iV ^ /i"7r+) 



F(M+ ^ ^i+N) 



TiN" e+TT") 



N 



TT^N) X 



F(iV^Z-7r+) TiN^^l+n- 



N 



■ N 



(7) 
(8) 
(9) 

(10) 



valid in rrii + < m^v < rUr — m^. Studying in subsection III A purely leptonic r-decays shown in Eq. ([5|, 
we will need the decay rates summed over all the standard light neutrino and antineutrino in the final state. The 
corresponding formulas take the form 



F(t e e+e vv) « {1 + 6n) 



I 



F(r- 
-F(r- 



e VeN) 



T{N e+e-^vi) 

TV ^ 
F(iV^ ^ e+e-vi) 



(11) 



■ N 



F(t — > e e /i vv) 



T[t- 
Sn ■ F(r 



e VeN) h On ■ F(t 



e-i^rN') 



Fw 
F(iV 



e^VeN) 



(1 + ME 



T{N - 



-+T{t- 
e+e^vi) 



F(jV^^eV-p^) 
Fat 

F(7V= ^ 



(12) 



■ N 



> e-VrN") 

T{t- ^l-UrN") 



e+e n) 



■ N 



valid in 2me < tojv < "^^t ~ "m-e- Here Sm = 0,1 for Dirac and Majorana case of sterile neutrino N, respectively. 
Summation in (11) and (12) runs over I = e,/i,r. The partial decay rates F(r, M — > XN) and T(N — )■ Yl) and the 
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total decay rate of sterile neutrino F^r involved in Eqs. (|7|-([l2|) are specified in Appendix. Implicitly all the partial 
decay rates include the corresponding threshold step- functions. For further convenience we rewrite Eq.(ll|, (12) in 
the form 



F(r- 



e e^e vi/) w (1 + Sn) 



FV 



(ei/Af)^(eeiy^) 



N 



N\ 



\U, 



Fjv 

\^+meN\^) 



(|t/rw|^ + \U^N\^\UrNf + W+ 1) | C/eW | ' | C/r JV | ' + 



(13) 



F(7 



where 



e e /i viy) 



■ N 



+ FV 



N 



■ N 



-((1 + 5N)a2\UrN\'' + (ai + 2(1 + SN)a2)\U^,N?\UrN? + (14) 



+ (^Arai + (1 + (5Ar)/3a2)|C/eivn(7,Ar|2 + {S^a^ + (1 + <5jv)/3a2) | f/ejv H t/^jv P + ai|t/ejv|^ + 

+ {l + &N)a2\Uf,N\''), 



/3 



N 



4.65, 



p(ej/iV)p(e/ii/) 
At J- at 



-p(eiy7V)-p(e/iiy) 



-p(/j!^Ar)-p(eei/^) ' 



a2 



r r r 



N 



p(/^i/Af)p(eeiy^) ■ 
It i TV 



(15) 
(16) 



In Eqs. ([13|)-([16|) we used notations r^^'^^FV''"'' introduced in Eqs. (|A.l[)-( |A.ll[ ). 

As we already mentioned, in the resonant regions of the sterile neutrino mass m^r, specified in Eqs. ([7])-([T2|, the 
intermediate sterile neutrino N, produced in r and meson M decays (see Figjl]), propagates as a real particle and 
decays at certain distance from the production point. If this distance is larger than the size of the detector, the sterile 
neutrino escapes from it before decaying and the signature of r — > Imr, r — > eelvv or M — > ttU cannot be recognized. 
In this case in order to calculate the rate of t or meson decay within a detector one should multiply the theoretical 
expressions F in ([7|-([T2]) by the probability P/v of sterile neutrino decay within a detector of the size Lb- Within 
reasonable approximations it takes the form [25) 



Pn w 1 — exp{—L£iT 



N), 



(17) 



where F^r is the total decay rate of sterile neutrino calculated in (A. 17). 

Then, the rates F^ of r and meson decays within detector volume should be estimated according to 



F^ = F X P, 



N, 



(18) 



where F are decay rates given by Eqs. ([7[)-( |12[ ). In our numerical analysis we take for concreteness Ljj = 10m which is 
typical for this kind of experiments. In Fig. 121 we plotted P/v v.s. sterile neutrino mass m^r for several values of mixing 
matrix elements |C//jvP- For illustration of typical tendencies we assumed in this plot |C/eAfp = — \Utn\'^- We 

do not use this assumption in our analysis. As seen, Pjv becomes small for tojv < 100 MeV even for rather large 
values of |J7|f^. Thus, in this region of rriN the effect of finite size of detector, described by Pn, significantly affects 
the decay rates of the studied processes and should be taken into account. 



III. LIMITS ON STERILE NEUTRINO MIXING UrN 



In the literature there are various limits on the mixing parameters UaN (with a — fj,, r) extracted from direct and 
indirect experimental searches [34j for this particle, in a wide region of its mass. A recent summary of these limits, 
extracted from the corresponding experimental data, can be found in Ref. [25] . In the present paper we focus on the 
least constrained mixing parameter UtN- In Fig. [sjwe show the exclusion plots for |[/t-jvP existing in the literature 
[25j together with our exclusion curves derived in the present section. For derivation of these curves we will analyze 
sterile neutrino contribution to the decays listed in ([5|-([6|. 

As seen from Eqs. ([7])-([l2]) the decay rates of the processes ^ depend on all the three UaN (with a — e,fi,T) 
mixing matrix elements. In the literature it is common practice to adopt some ad hoc assumptions on their relative 
size in order to extract limits on them from the experimental bounds on the corresponding decay rates. In particular, 
limits from CHARM [35] and NOMAD :363 plotted in Fig. [s] assume It/rAfl 3> |f/eAf|, |C^/iAr|- These assumptions may 
reduce reliability of the obtained limits. Below we derive analytic expressions for limits on iC/rAfP in different mass 
ranges of toat without any kind of such assumptions. 
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FIG. 2: The probability Pjv of sterile neutrino decay within a detector of the size of Lo = 10 m versus sterile neutrino mass 
mjv for several values of mixing matrix elements |i7ijvp, assuming |f/eivp = If^M^vl^ — If/rjvl^. 



A. Purely leptonic decays 



First we exploit for extraction of | L/t-w | the following experimental results for the branching ratios of purely leptonic 
T-decays [31] 



Br(T — )■ e e^e v^Vr) 
Br{T^ — fi^ Vfj^Vr) 



^ (2.8 ± 1.5) X 10"^ 
< 3.6 X 10"^ 



(19) 
(20) 



The first decay has been observed experimentally and its experimentally measured branching ratio agrees with the SM 
prediction within the standard deviation A'^^^'(r" e~ vv) = 1.5 x 10~^. Neutrino assignment in the final states 



of the decays (19l-(20) corresponds to what is suggested by the SM. However, in the experiments, measuring these 



decays, the final state neutrinos cannot be actually identified. Therefore, considering beyond the SM mechanisms 
with LFV one should take into account the possibility that all the light neutrinos v^, Vfi and Vj- may contribute to 



the final state of the decays (19)-(20). Formulas p3|-Jl4h were derived for the very this case. They describe the 



sterile neutrino resonant contribution (diagram FigTTra)yto the decays (19)-(20l and will be used in the analysis of 
this subsection. 

We also assume that the sterile neutrino contribution to the process (19), if exists, should be less than A'^^p. For 
the decay (20), not yet observed experimentally, there exists only the above indicated upper bound and the sterile 



neutrino contribution has to obey this bound. 

Taking into account the finite detector size effect according to Eq. ( 18 1 we write for decay rate F^ within detector 
volume 



r^(T- 



e e^l vv) 



r(7 



'I vv) X Pj 



(21) 



with r(T^ — > e^e^l^vv) given by (11), (12). As we discussed in the previous section, the probability P/v of sterile 
neutrino decay within detector becomes rather small for mjv < 100 MeV. Therefore, in this mass range we may 
approximate the expression in (17) by ~ LjjTn. This is a reasonable approximation for this part of our analysis 
since the limits, which will be obtained here, correspond to the exclusion curve (a) in Figs. [3] and curves in Figs. |4j [5] 
located in the region rriN < 100 MeV, where LdTm ^ 0.01. 
In this approximation we find from (13) and (21 1 



i AT ^ 



D' 



\UrNf + \U^Nf\UrNf + (/? + iWeN^PrN^ + (22) 



According to our assumption, discussed after Eqs. (19)-(20), we require 

r^r(r^ vv 



< A^^P(r" ^ e-e+e-'w) « 1.5 x 10" 



(23) 
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where = (290.6 ± 1.0) x IQ-^^ s is the r-lepton mean hfe |34) . Then we obtain the following upper limits 



A'^^p(t — > e e+e vv 



y (1 + on)Ld Tr 



\UrNU^N\ < 



|f/rwC/eJv| < 



A'=^p(t" 



e e+e i^u) 



(/3 + 1) r-^r^^''-(l + ,5Ar)iz3 



(25) 



Similarly, we derive limits based on the experimental bound (20). Using Eq. (14), we find 




Br'^^P{T — !■ fi e+e lyiy) 



/i e^e z^j/j 



(TV' 



N 



Br''^P{T- 



/i e+e vu) 



(e!^Af)-p(e/jj^) 



+ r^'^-^'ri^"'^-' )(<5Arai + (1 + <5w)/3a2)Lc ' 



(26) 
(27) 
(28) 



Here, Br'^^P denotes left-hand side of the experimental bound in (20). The limits (24)-(28) are plotted in Fig. 3][^ 
for the case of sterile Majorana neutrino. Drawing the exclusion curves, we selected the most stringent limit among 
(24)-([28| for each mass value rriN within the studied mass range. As seen, the present experimental data (19)-(20) 
on purely leptonic r-decays set rather weak constraints on |?7rAr| and on IUtnUcnI, lUrNUf^Nl in the mass region 
IMeV < mj\i < 100 MeV. Our limits on |[/rjv|, corresponding to the curve (a) in Figlsj are significantly weaker than 

TNUeNlt 



the limitations from other searches shown in Fig. [sj However, our limits for \Ur 
our best knowledge are new in this mass region. 



and \UrNU^N\ in Figs, mmto 



B. Leptonic and semileptonic decays 



Now we combine the purely leptonic r-decays considered in the previous subsection with the semileptonic decays 
of r and K, Z?- mesons using the experimental data ( 19 ) , ( 20 ) and the experimental limits on the following branching 
ratios 1341: 



Br{T- 
Br{K+ 
Br{D+ 



-) < 1.2 X 10"^ 



7r"e+e+) < 6.4 x 10"^°, 
7r~e+e+) < 3.6 x 10"^ 



Brir- 7r"7r+/i+) < 7 x 10"®, 
Br{K+ TT+fi-e+) < 1.3 x 10"^^ 
Br{D+ n+fi-e+) < 3.4 x 10"^ 



(29) 
(30) 
(31) 



Assuming that in all these decays sterile neutrino N contributes resonantly wc should limit ourselves to the mass 
domain: 



245 MeV < mjv < 



1637 MeV. 



(32) 



Within this mass domain the experimental bounds ( 30 ) contribute to our analysis only up to 
TOjv < rriii — TOe ~ 493.2 MeV corresponding to the mass range of the resonant contribution of sterile neutrino 



to these decays of K-meson. In the above list (29)-(31) one could also include the existing experimental bounds 



on the other LNV and LFV decays of r and Ds,B mesons. However, they have negligible impact on our results 
presented below. 

In this part of our analysis we put Pn 



1 for the probability (see Eq. ( 17 )) of decay of nearly on- mass-shell sterile 



neutrino, resonantly contributing to the analyzed processes. Thus we assume that these processes occur completely 
within a detector volume. This is a good approximation for the case of the limits on UrN^ which will be derived here 
and displayed in Fig. [3] as curve (b) . To see this one can check the plot for P/v shown in Fig. |2] 

In the mass domain ([32]) we can use Eqs. ([7l)-([l4| for the corresponding decay rates. Below we combine these 
formulas in a system of equations. Solving them with respect to \UrN\ and applying the experimental bounds (19), 



(20) and (29)- (31) we find upper limits on this mixing parameter. For our purpose it is sufficient to use ether of the 



two experimental bounds (19), (20). We select (19) which leads to a bit more stringent limits on IC/tAtI- 
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Let us introduce the following notations 



Fee{T) 
FeeiM) 



■ On) i- T i ^ 



(1 



Br''^P{M+ 7r~e+ 

p(eAf)p(e7r) 



(7rJV)p(i7r) ^ 



(33) 



M 



N 



(r 



(eW)-p(p7r) 



M 



N 



) TAf 



where t^,tm are mean lives of r and M = K'^,D^; the right-hand sides of the experimental bounds in (29)-(31) are 
denoted by Br^^P; the quantity A'^^p was introduced after Eqs. (19) and (20). 



Now we can rewrite the experimental limits on (19) and (29)-(31) in the form 



ae\UeN\^ + a^lf^MJVp + ar\UrN\ 



ae\UeN\^ + a^\Ufj.N\^ + ai-|[/rjvp 



r < F^i{t), 
<F,AM). 



Here I = e,^. Solving (34)-(36) we find 



\UrN? < CiF^eir) + C2F^e{T) + C3F^^(r) + C4Fee(M) + c^F^^iM). 



where 



Ci = Ot, C2 = flg — 2ar! C3 = flp 



C4 



(/3 - l)ae - /3a^, C5 = (/3 - 1)0^ 



(34) 
(35) 
(36) 

(37) 
(38) 



We have checked that in the mass region (|32| all the coefficients Ci > 0. The parameter /3 is defined in (15). We 
plotted the corresponding exclusion curve in Figjs] labeled by (b) for the case of Majorana sterile neutrino. As seen, 
our limits are more stringent than the existing ones from CHARM }35j and DELPHI [37. experiments in the sterile 
neutrino mass region 300 MeV< toa? <900 MeV. Note that in difference from the existing limits on \UrN\ our limits 
are model independent in the sense that we have not made any assumptions on the other two mixing parameters 
\UeN\ and |C/mjv|. Instead, we excluded them combining the experimental limits on the branching ratios of different 



processes (19), (20) and (29 )-( 31). 



Sterile neutrino in the final state 



Other experimental data which we apply for deriving limits on UrN are |34) 

Br{T- iDiUr) = (17.85[17.36] ± 0.05)%, 
Br{T-^TT-Ur) = (10.91 ±0.07)%, 



(39) 
(40) 



where in the first line the central value 17.85 corresponds to ^ = e and 17.36 to I — fi. Both these experimental results 
agree with the SM predictions within the standard deviations A^^p{t Ivv) = 0.05% and A'^^p(t — lirv) = 0.07%. 
We already commented in subsection IIIA (after Eqs. (19), (20)), that in the reported experimental results like in 
Eqs. ([l9|)-(20) and (39)-(40) the final state neutrino assignment t^e,^,r is made according to what is suggested by the 



SM. However, in the experiments, measuring these decays, the final state neutrinos cannot be actually identified and 
are observed as a missing energy signature. Therefore, it is liable to imagine that instead of one or even both of the 
standard light neutrinos in the final states of decays in (39)-(40) there may occur some other neutral particles such 
as sterile neutrinos. We assume that in these modes of r-decay appears one sterile neutrino N accompanied by any 
of v^tV^^Vt. Its mass must satisfy to mjv < — mi and uin < rrir — ra^^ for the decays (39) and (40) respectively. 
We also assume that this contribution, if exists, should be less than the corresponding stanHard deviation A'^^p. 



The contribution of sterile neutrino to ( 39 ) , ( 40 ) in the form 



-N 



(41) 



should be less than the corresponding A'^^^ since (39), (40) are in agreement with the SM. 
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Therefore, using (A. 2 1 and (A.3I we find the hmits 



\UrN? < Min J ^ 



p(TrJV) 



(42) 



where the minimal of the two values in the curl brackets are selected for each value oiniM- The corresponding exclusion 
curve is shown in Fig 3] and comprises the two parts (c) and (e). The part (c) is dominated by the constraints on purely 
leptonic r-decay mode while the part (e) is mainly due to the semileptonic mode shown in (41 ). The exclusion curve 
(c), (e) cover a mass region < mjy < to^. — TO^r w 1640 MeV. This curve sets new limits on Utn for < m^r < 70 
MeV and 300 MeV < < 700 MeV. In the region 500 MeV < niN < 700 MeV they are less stringent than our 
limits derived in the previous subsection from the data (19), (201, (29)-(31| and corresponding to the curve (b) in 
Fig. |3] For m^r < 100 MeV the part (c) of our exclusion curve is nearly constant and our limits for this mass range 
can be displayed as 



\Ur 



N\ 



< 2.9 X 10" 



for < < lOOMcV. 



(43) 



As we discussed previously, sterile neutrino produced in (41) can decay within a detector with a probability Pi 



N 



defined in (17). This would result in appearance of a displaced vertex attributed to this sort of decay in addition 



to the production vertex (41). The limit in (42) does not take into account such a possibility and sum up the event 
rates of sterile neutrino decay both within and outside a detector. However, one can imagine an experiment where the 
displaced vertices of the above mentioned type are looked for and are either observed or, more probably, excluded at 
certain confidence level. For the latter case our limits in the region rnjv >100 MeV would drastically change. In order 



to illustrate the influence of this additional criterium of event selection on our limits we impose on the processes (41 ) 
a condition that sterile neutrino decays outside detector. This results in multiplication of the corresponding decay 
rate formulas (A. 2), (A. 3) by the probability factor 1 — Pm- The modified limits take the form 



IC/rATp < Min^ ^ 



vvl 



(ttN) 



X cxp{LdT%). 



(44) 



Here we us ed an i nequal ity exp (LdTn) < exp (LuT^), where F' 



defined in (A.17) 



N 



(A. 18). In this case our exclusion curve for \UrN 
(|42|) changes its part (e) to (d) leaving the part (c) intact. Now the exclusion curve (c)-(d) covers a mass region 



- ae{mN) + a^{mN) + ariniN) with ae,^,r 
^ in Fig. 3 in comparison to the case of 



< niN < — TOtt ~ 180 MeV. Note again that this is just an illustration of an impact of as yet non-existing ex- 
perimental data allowing discrimination of the events with the displaced vertices associated with the sterile neutrino 
decay. 



IV. SUMMARY AND CONCLUSIONS 



We studied resonant contribution of sterile neutrino to leptonic and semileptonic decays of r as well as to some 
semileptonic decays of K and D mesons. Comparison of our predictions with the corresponding experimental data 
on these decays allowed us to extract new limits on the mixing matrix element Utn shown in Fig. |3] as curves (b) , 
(c), (e). In the two domains of the sterile neutrino mass < mAr < 70 MeV and 300 MeV < < 900 MeV our 
limits on UtN are more stringent than the limits existing in the literature. For < ttiat < 100 MeV our limit to a 
good approximation is |C/rArP £ 2.9 x 10"'^. We also obtained new, although not stringent, limits on the products 
\UrNUeN\ and \UrNUfj_N\ showu in Figs. |4][5] To our best knowledge there do not exist in the literature the limits on 
these products of the mixing matrix elements for ttin < 100 MeV. 

Our limits derived from the experimental results (39)-(40) are, to certain extent, conservative estimates. In fact, let 
us assume that in derivation of these experimental values were applied specific kinematical criteria for event selection, 
suppressing possible contribution of t — ^ vriV, liyN-decays with a massive neutral particle N, such as sterile neutrino, 
instead of the nearly massless neutrino. Then, taking properly into account these criteria in derivation of limits on 
|{7rAr| would have to strengthen them in comparison with our limits in Fig. |3] This sort of analysis is beyond the 
scope of the present paper and requires many additional and imknown for us details on the derivation of (39l-(40) 
carried out by the corresponding experimental groups. 

We consider as an important point of our analysis its model independent character in the sense that we do not refer 
to any sort of ad hoc assumptions about other two mixing matrix elements UeN and U^n- Such assumptions are typical 
for the existing literature on this subject. In particular, the limits of CHARM [35] and NOMAD ,36J collaborations 
shown in Fig. |3] were obtained under the assumption |t7T-jv| 3> |C^/jAf|, It^ewl- At first site this assumption looks 
reasonable since the existing limits on |t/;^Ar| and \UeN\ are very stringent (see, for instance, Ref. |25|). However, they 
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FIG. 3: Exclusion curves for |(7rjvp from the present analysis, denoted by (a)-(e), and the exclusion curves existing in the 
literature derived from CHARM '35', NOMAD [36 and DELPHI [S^ searches for sterile neutrino decays. The latter curves 
are taken from Ref. ^25) . 
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FIG. 5: Exclusion curves for \UtnU^n\ from the present analysis. 

were also obtained under the assumptions of this type. To our mind these observations should be taken into account 
in assessment of the limits on the sterile neutrino mixing matrix elements UaN- In some cases these limits may be 
rather stringent mainly because of this sort of assumptions. 
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Appendix: Partial decay rates 

Here we specify the partial decay rates involved in Eqs. ([7l)-(|l2|). For more details and discussion we refer reader 
toRefs. [211 [21 Hi] 

The decay rates of mesons and r to the final states with sterile neutrino N: 

T{M+^l+N) ^ \U,^\^^fl^\VM\'mljXHxlx%A){xj+x%-{x',-x%f)^ |f/.Arpr(^^\ (A.l) 

on 

TiT-^TT-N) = \UrN\^^mlf^\Vud\'Fp{zN,zp) = \UrN\^ri-''\ (A.2) 

Tir-^l-i^iN) - \U,^\^-^mlh{zN,z,„zi) = \UrNM''''\ (A.3) 

TiT-^l-,.rN) = \UiN\^-^mlhizM,z,^,zi) = \UiNM''''\ (A.4) 
Here we denoted z,; — nii/m^, Xj = mi/niM with to^ = m]^^mp,mi. The kinematical functions Fp(x,y), Ii(x,y,z) 



are defined in (A. 13) 
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The partial decay rates heavy sterile neutrino, N including leptonic and semileptonic decay modes. In the latter 
case the final hadronic states for low neutrino masses rriN < rnp is represented by the lightest mesons while for larger 
rufq > nip by gg-pairs as suggested by Bloom-Gilman duality [35] • This inclusive approach [3T] allows one to reduce 
uncertainties in the leptonic decay constants /m of mesons starting from p-meson, some of which are only known in 
phenomenological models (for more details see [21] )• The list of the sterile neutrino decay rates is as follows: 



r{N^i^i+^i,) - 

r{N^,yiJ^l+) = 



r(iv^zj-p+) 

V{N l^ud) 
TiN ^ ivi, qq) 



2 ^5 _ irr |2p(3l^) 



647r 



G 



(A.5) 
(A.6) 

(A.7) 

(A.8) 
(A.9) 
(A.IO) 
(A.11) 



Here P — tt,K. The decay constants are /^r = ISOMeV, Jk = 159MeV. We denoted yi — mi/niM with rrii ~ 
mi,mp,mq. The CKM factors in Eq. (A.8 1 is Vjr = V^J^^^ , Vk = V^J^'^^ . For the quark masses we use the values 
m„ « m<i = 3.5 MeV, = 105 MeV, rric = 1.27 GeV, rub = 4.2 GeV. In Eqs. ( |A.10[ ), ( jA.lip we denoted u^u,c, t; 
d = d,s,b and q — u,d, c, s, b, t. The SM neutral current couplings of leptons and quarks are 

g[ = -l/2 + sm^9w, gl = 1/2 - (2/3) sin^ 0w , gt = -1/2 + (1/3) sm^ 9w, 

gl = -(2/3) sin2 Ow. g'k = (1/3) sin^ dw 



(A.12) 



gp = siii^ew, 



The kinematical functions in Eqs. (A.l)-(A.ll) are 



ds 



h{x,y,z) = 12 / -{s-x^-y'){l + z^-s)X'/\s,x^y^)X'/\l,s,z^), 



ds 



h{x,y,z)^2Ayz / —{1 + - s)\^'^{s,y\z')\^/\l,s,x^), 



(A.13) 



(A.14) 



Fp{x, y) - Ai/2(i, x\y^)[{l + x^){l + - y^) - Ax% (A.15) 

The total decay rate Fjv of the heavy neutrino TV is equal to the sum of the partial decay rates in Eqs. ( A.5)-( A.11[ ), 
which we write in the form: 



Tn = ^ (1 + 5n) [r(7V ^ l^n+) + T{N ^ i^itn,)- 



(A.16) 



ii.h.n 



where we denoted the hadronic states V.'^ — ,du,su,dc,sc and = P^,qq. We introduced the factor Sn = 1 
for Majorana and 6^ = for Dirac neutrino N. Its appearance is related with the fact that for Majorana neutrinos 
both charge conjugate final states are allowed: N — >■ l^ I2 I'l, , l^ I2 Pi,; N — > Z^Z^i^;^ , ^^/^i/;^ and N — > PV,^ . For 
convenience we write Eq. (A. 16) in the form: 



Fat = aeirUN) ■ \UeN\^ + CLfiimN) ■ \U^n\'^ + arimN) ■ PrNl"^ 



(A.17) 
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where 



ai{mN) = (1 + 5n) 



with I, I2 = e, /i, T. In the inclusive approach the hadronic contribution is calculated as 



(A.18) 



(A.19) 



The parameter /ig denotes the mass threshold from which we start taking into account hadronic contributions via qq 
production. In Refs. [HI [52] we have shown that the reasonable choice is = mp+ = 775.8 MeV, which we also use 
in the analysis of present paper. In Fig|6]we plotted r^vo = ^NiUeN — U^n — UtN = 1) as a function of the sterile 
neutrino mass toat. 
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FIG. 6: Sterile neutrino decay rate Fjv for the particular case of Uen = f^piv = Utn ~ 1 denoted by Fivo- 
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